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Abstract 

We investigate the effect of a large number of kinks on the gravitational power radiated by cosmic 
string loops. We show that the total power radiated by a loop with N left-moving and right-moving kinks 
is proportional to N and increases with the typical kink angle. We then apply these results to loops 
containing junctions which give rise to a proliferation of the number of sharp kinks. We show that the 
time of gravitational decay of these loops is smaller than previously assumed. In light of this we revisit the 
04 ■ gravitational wave burst predictions from a network containing such loops. We find there is no parameter 

5^ t regime in which the rate of individual kink bursts is enhanced with respect to standard networks. By 

contrast, there remains a region of parameter space for which the kink- kink bursts dominate the stochastic 
background. Finally, we discuss the order of magnitude of the typical number of sharp kinks resulting from 
kink proliferation on loops with junctions. 
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1 Introduction 



^ ' Cosmic strings arc one-dimensional topological defects of cosmic size originally introduced by Kibble pQ which 
• may have been produced during spontaneous symmetry-breaking phase transitions in the early Universe. Rc- 
\ cently, it has been realized that fundamental objects from string theory could, in certain scenarios, expand to 
cosmic size and also play the role of cosmic (super-)strings (see [2] and references therein). A very promising 
strategy to detect them relies on their gravitational wave emission, and in particular on their emission of high 
t> ■ frequency bursts at cusps and kinks [31 H] which could be accessible to the advanced versions of the network 
\ of ground based interferometers LIGO- VIRGO and to the future space interferometer LISA in a large region 
of parameter space. Although cosmic strings have never been observed so far, a lot of effort is made in this 
| direction since their detection would provide unique insight on either high energy physics or string theory. 

Several authors have studied how certain properties of cosmic supcrstrings could help distinguishing their 
gravitational wave burst signal from that of field theory cosmic strings, such as their reduced reconncction 
probability [5], their reconnection process itself [5] or the fact that they can move in extra dimensions [7]. 
Another feature of supcrstrings that drastically affects their gravitational wave burst signal is the existence of 
junctions between strings. A cosmic superstring network contains F- and D-strings as well as bound states made 
up of both types called (p,q) strings. The crossing of two strings of different types can lead to the creation of 
a new segment of string joining the original ones at two vertices also called Y-junctions [SJ O [TUl HU E2J. Note 
that junctions can also be formed in specific models of symmetry-breaking [131 HH [T5] . 

The impact of junctions on the gravitational wave burst emission from a cosmic string loop was investigated 
in [THl [T7] . It was shown that the most important effect they induce is the exponential proliferation of large 
amplitude kinks. In [TS], the effect of kink proliferation on the signal emitted by a whole network of loops 
was investigated. It was found that, for models where the fraction q of loops that contain junctions and 
the average number k' of sharp kinks on these loops satisfy the condition qk! > 1, then the observational 
predictions differ from the standard case examined in [31 119j , mainly due to the presence of a strong stochastic 
background resulting from the incoherent superposition of bursts produced at kink-kink collisions potentially 
hiding individual bursts in certain regions of parameter space. 

One of the key assumptions in this analysis was that the lifetime of the loops with junctions (or of their 
daughter loops in the case where the junctions unzip) was unaffected by the large number of kinks on them. To 
be more precise, it was assumed that the gravitational power radiated by these loops was given by the standard 
formula P = TGfi 2 with T ~ 50, a result first obtained by Vachaspati and Vilenkin in [201 and a priori only 
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valid for rather smooth loops. This rate of energy loss determines the lifetime of the loops in the network and 
therefore enters crucially in the estimation of the density of loops at any time and ultilmately in the rate of 
observed events. 

In this paper, we investigate the validity of this assumption by computing the radiation from piecewise linear 
loops with a large number of sharp kinks. The radiation from kinky loops was already explored in the past, 
but always with a different focus. The first calculations date back to the work of Vachaspati and Garfmklc 
(21] who showed that the power radiated by kinky cuspless loops was finite in every direction and computed 
analytically the power radiated by a particular class of loops with 4 kinks in total. Other simple configurations 
were studied in [22] by Allen, Casper and Ottewill. In [23], Allen and Casper derived an analytical expression 
for the total power radiated by any piecewise linear loop made of N straight segments as a sum of 0(N 4 ) 
polynomial and logarithmic terms. Their idea was to provide a way of computing the power from smooth loops 
by approximating them by increasingly complex polygons with angles at the vertices tending to zero, a situation 
completely different than the one we arc considering here since we are interested in a large number of sharp 
kinks. Although their formula can be used to compute the power from any given loop with many sharp kinks, 
its complicated form does not allow to determine analytically the generic large N behavior for such loops. 

Here, we develop a completely different approach, based on an analytical averaging of the power over loop 
trajectories with N sharp kinks leading to an expression where the large N behavior can be extracted. This is 
the first time that the power radiated by loops with many sharp kinks is investigated using a fully relativistic 
formalism, namely Weinberg's formula and our main result is that the power radiated is proportional to the 
number of kinkfl 

We then apply our results to loops with junctions and revisit the estimate for their number density and the 
gravitational wave signal predictions made in |18] , In particular, we investigate the modifications to the rate 
of individual kink bursts and the amplitude of the stochastic backgrounds formed by the superposition of kink 
bursts and of kink-kink bursts. 

Finally, our results provide a new step towards understanding how long kink proliferation can last and 
towards estimating the resulting value of the parameter k' . Indeed, the expression of the power radiated as 
a function of the number of kinks enables us to compute the energy lost as kinks proliferate on the loop and 
therefore to estimate the time it takes for a loop to decay by gravitational radiation. We also discuss qualitatively 
how other mechanisms such as intercommutations, unzipping of junctions and gravitational backreaction affect 
our estimate. 

The paper is organized as follows. In section [21 we study the radiation from piecewise loops with a large 
number of kinks. In section 3, we use our results to revisit the observational predictions taking into account 
the shorter lifetime of the loops with junctions. In Section 4, we discuss the implications of our results on the 
duration of kink proliferation and therefore on the number of kinks. Section 5 contains our conclusions. 

2 Radiation by a loop with many kinks 

In this section, we compute the power radiated from periodic piecewise loops with N kinks. Sec. 12. H is devoted 
to the definition of the class of loops that we consider. Then, in Sec. 12.21 we use Weinberg's formula to derive a 
general expression for the power radiated by the loops in this class. In Sec. 12. 31 we analytically investigate the 
dependence of the power on the number of kinks. This section contains our main result, namely the fact that 
for large N, the power becomes proportional to N. In Sec. 12.41 we use numerical calculations to investigate 
how our result depends on the amplitude of the kinks. Finally, in Sec. 12.5] we comment on the applicability of 
our results to the case of loops with junctions. 

2.1 Piecewise constant loops: the setup 

Let us start by setting the notations and defining the class of loops that we use. We consider a cosmic string loop 
whose dynamics is governed by the Nambu action. We assume flat spacctime background geometry ( — h H — h) 
(we only consider loops small compared to the horizon size) and we use the standard conformal time gauge 
to parametrize the string's worldsheet so that the string is described by its spatial coordinates x(tr, t) where t 
coincides with Lorentz time. The gauge constraints can then be written (with ' and ' standing for derivatives 

1 We note that the linear dependence of the power on the number of kinks had already been proposed in [24] but the argument 
was based on the quadrupole radiation formula which only applies to slowly moving sources and whose validity is therefore not 
clear in the context of cosmic strings. 
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with respect to a and t) 



x' • x = 0, 



± z = 1. 



The wave-like equation of motion leads to the form 

1 



x(tr, i) = - (x+ (<r + t) + x (a - t)) 



(1) 
(2) 



(3) 



with = xj 2 = 1 in order to satisfy the gauge constraints (here ' denotes a derivative to the only variable of 
x + or x ) . 

As can easily be checked, the total power radiated will be independent of the invariant length of the string L 
so, in order to simplify the expressions, we set L = 2tt so that at any time, a G [0, 2tt]. The evolution is periodic 
and x + and x + are periodic functions 



x±(0) 



x±(2tt). 



(4) 



The definition of these two functions on the interval [0, 27r] is set by the inital conditions (x'(t = 0,tr) and 
x(t = 0, a)). The whole evolution of the loop can then be determined trivially using periodicity. We now define 
the class of initial conditions that we use. 



Let A± jT1 , n = 0..N — 1 be two sets of N unit vectors such that ^2n=o A±. n = 0. We can define the derivative 
of the functions x±(u) as a piecewise constant functions: 



x' ± (it) = A± >n for u G 



2nn 2(n + l)n 



(5) 



By integration, x± arc picccwisc linear functions given by (the integration constants vanish in the center of 
mass frame of the loop) 



x± (it) 



2tt 
A 



2nn\ 



for u G 



2mr 2(n+l)7r 



N 



(6) 



(with the convention X) 9 =o A±. g = 0). The condition (J3| is clearly satisfied since ^2^=0 -A-±, n = 0. 
A loop defined like this has ./V left moving and A right moving kinks evenly spaced on the string^ 



2.2 Radiated power 

Because the evolution is periodic the frequencies of emission are all multiples of the fundamental mode of the 
string ui m = mwo where ujq = = 2. The total power radiated by the loop per unit solid angle in the direction 
given by the unit vector k is a sum over all such frequencies and can be computed using Weinberg's formula 

m 

HP °° Ci) 2 m 2 

5r = E ^^^. lv {V)T^{V m ^ m )T^{\, m ^ m ) (7) 

OK L — ' 7T 

m— 1 

where k m = oj m k, r* J '(k, w) is the spacetime Fourier transform of the stress-energy tensor of the string and 
Ajj ) i p (k) is a projection operator defined as Ay.; p (k) = PuPj P — \PijPi m where Pjj = 5^ — kikj. Note that 
Py /c l = 0: Pij is the projector on the plane perpendicular to k. Here S is the Kronecker symbol and all indexes 
i,j,l,p range from 1 to 3. 

The (spacetime Fourier transform of the) stress-energy tensor can be conveniently written under the factor- 
ized form [3] 

T«(k m ,w m ) = Jlj£(k,m)i! ) (k,m) ) (8) 

2 We are mainly interested in such symmetric configurations because we expect any mechanism responsible for the creation of 
many kinks to be symmetric in the left and the right movers. In particular, kink proliferation on loops with junctions that we 
discuss in Sec. 12.51 creates similar distributions of left-moving and right-moving kinks. We will however comment on the power 
radiated by loops with N 3> 1 left-moving kinks and only a small number of right moving kinks at the end of Sec. 12.31 
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where ^ denotes symmetrization with respect to the indexes i and j and 



p Z7i 

F±(k,m)= x'l{u)e^ muia{±u -^ ±{u)) du. (9) 
Jo 

The integrals above can be easily computed using and and yield 



P± = A j2 ci )n e^( ±n - k -(s::o A ±.0) (io) 



n=0 

where we defined 



.4 



k • A±,„ T 

Both I± and the C±. n depend on k and m but from now on, in an effort to avoid heavy notation as much 
as possible, we leave this dependence implicit. Note that „ is not singular at k = ±A± y1l , and there- 
fore I± is not singular at any k = ±A± -n . More precisely, if k = A +no , we have /+(k = A + ,„ ,m) = 

^-A i + no e 2l ^ 1 ( no '' U '^-=o A+ s )) + O(i). Therefore, although I + is finite in those directions, it doesn't tend 
to zero as m goes to infinity. However, the series in ([7]) still converges because the operator A projects out this 
first term since it is proportional to A +i „ = k. As a consequence, 4? is finite in every direction in perfect 
agreement with the findings of [3T] in the case N = 2 (except in the degenarate case of a persistant cusp when 
one of the A +jra is exactly equal to one of the A_, 9 ). 

Finally, the total power radiated by the loop is the integral of (J7]) over all directions 

(12) 

2.3 Dependence on the number of kinks: analytical study 

2.3.1 Loop distribution 

We now would like to average over the class of loop trajectories defined in Section 12.11 for fixed values of the 
number of kinks iV and study how the total power radiated varies with N. In order to do so, we first need 
to specify the weight of each loop of the class or equivalently, a random procedure to construct such loops. 
This raises the question of the naturalness of the probability distribution that we will use: what do realistic 
kinky loops look like? Are the kink angles typically large or small? For now, we do not address these complex 
questions (we will focus starting from Sec. 12.51 on a particular mechanism creating many kinks, namely kink 
proliferation on loops with junctions) and make the arbitrary choice of using a "flat" distribution (in a sense 
that we explain in the next paragraph) for which the average can be computed analytically. In Sec. 12.41 we use 
numerical computations to discuss how different distributions can affect the result. 

The simplest idea to obtain a flat distribution is to draw randomly each of the A + ,„ and A_ ; „ for n = 0..N—1 
with uniform probability on the unit sphere. However, we must satisfy the conditions ^2^=0 -^-±-« = 0- ^ n an 
attempt to generalize the trajectories used in [2T] in the case N = 2, we use the following construction: we 
restrict our attention to even values of N and draw each of the A +i „, n — Q..N/2— 1 and A_. n , n — 0.-^/2 — 1 with 
uniform probability on the unit sphere. For n = N/2..N — 1, we then define the remaining A±, n = — A± , n -N/2 
so that the sums vanish. To be precise, this yields a flat distribution but only over a subset of our initial class. 
Though such a construction is probably the simplest, other are possible and the robustness of our results to the 
distribution of loops will be discussed in the next section. 

2.3.2 Average power 

Since there is no preferred direction in our loop distribution, we expect to be independent of k. We can 
therefore set k = (1, 0, 0) from now on to simplify our calculations. From J7J) and (|BJ), we have 



dP\ Gfi 2 ±5° 



rfk / 47T 3 

m— 1 



£ »rA„„((/y/'' ■ /;•/") (l l + l*_ + illL)) (13) 
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where from now on the brackets (■) denote an average over the class of loops. Since the A + and A_ are drawn 
independently, the averages over I + and /_ factorize and we are left with 

% ) = e ™ 2a «* (<4- 4> {p-i p ~) + (i 3 -i 1 -) + (4*4) + (4*4) 

' m— 1 

(14) 



Therefore, the building blocks of our calculation will be (I^I^ J and (^P_ P* j which, as we explain now, are com- 
plex conjugate quantities. Indeed, by remarking that (with obvious notations) C!_ n (A_. n ) = (C\_ n (— A_. n ))* 
and P_ ({A_.„}) = — Jl ({— A_ -n })* we can immediately conclude that 



1-1°-) = (I- ({A_. n })/I({A_, n })*) (15) 

4({-A_,„})*4({-A_,„})) (16) 

7;({-A_, n })4({-A_, n })*)* (17) 

44*\* (18) 



where in going from the third to the last line we used the fact that any set of vectors {A_ irl } and the set of 
their opposites {— A_ jrl } has the same probability with our procedure. 



We will now focus on one of these objects, say and drop all the indices + so that we now use the 

notation instead of C\ n and A„ instead of A + .„. Because the calculation is somewhat long and technical, 
we save all the details for Appendix E] and simply present the outline of the calculation here before jumping 
to the result. Starting from formula ([TO]), it is clear that (l+4*) will be a sum of average values of products 

of terms of the form C^C^I with several exponentials of scalar products of the type k ■ A s . The first step is 
to factorize these average values as much as possible using the fact that the random unit vectors A n for n 
between and N/2 — 1 are drawn independently to obtain a product of expectation values each involving only 
one of the vectors A„. We can then use the symetries to reduce the number of components that we need to 
calculate. Indeed, since we had set k = (1, 0, 0) and we will eventually contract (44*) with the operator A^jp 

which projects away any term proportional to k, any component of (44*) bearing an index 1 will disappear. 
Furthermore, after remarking that directions 2 and 3 are equivalent, it becomes clear that the only relevant 
quantity to compute is (44*) = (44*) which can he expressed in terms of the three following functions 
computed in Appendix [B] 

e — kA ") = - gin — — , (19) 



27T771 V N 



M - (clC^ N/2 e^^) = cos - _ sin , (20) 

and finally 

/^2^2*\ n. /4?rm\ N . /4?rm\ 
a = (0*0?) = Cm ^_j + _ Bm ^_j_l | (21) 



where the special cosine integral is defined as 



Cin(z) = f X 1 — dt. (22) 



Given our procedure to draw the vectors A„, it is clear that none of the average values defined here actually 
depends on the index n. Also note that all three quantities a, [i and s are real. 
Now jumping to the result of the calculation, we have 

(44*) = ^2 ( a + (-l) Tn s W/2 ~V) - S a ) + Q» (23) 
where the tensor Q lJ will be projected away by Aiji p . In other words, we can now replace in our calculations 
44*) ^ ^ e di a g° na l tensor 

N 

7T1 2 










) i 


1 


• 










(a + (-lr^-V) 1 . (24) 
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ri TV 



Note in particular that because the part of \I\Pj* ) that is not projected away is real, we can in practice also 



replace j in (jM|) by the expression given in 

When performing the contraction in (|14|) , the only term that appears is therefore the square of the factor in 
flUD and it appears 8 times. This leaves us with 

(%) = ¥ 2 £ [a(m,iV) + (-l)™ S (™,7V)^-V(m,7V)] 2 (25) 

* ' m— 1 

where we now indicated explicitely the dependence on m and TV. Note that all three quantities <r, fi and s are 
only functions of ^ so from now on we will slightly abuse notations and write cr(ra/TV) for example. 

2.3.3 Large TV behavior 

Our goal is now to understand the large TV behavior of (|25p , for which it is convenient to expand the square of 
the brackets, thus expressing (4^) as the sum of three terms 

dP\ 2Gu 2 f±^N 2 o fm\ TV 2 ,/m\ /m\N-2 _ +2? _ TV 2 /m\ /m\ /toN^/ 2 " 1 ^ 




(26) 

We will now show that the dominant term in the TV — > oo limit is the first one. The behavior of this 
term is simple to obtain because the dependence on m and TV is entirely of the form m/N so we can use a 
Ricmann integral to show that it is proportional to TV. For the other terms, precisely determining an equivalent 
is much more complicated but we just need to show that they are subdominant. For the sake of simplicity, 
we chose to present here only qualitative arguments for this instead of a complete but much more technical proof. 



Let us start with the first sum. In the large TV limit, we clearly have 

+ °° at-2 . . p oo 

E§- 2 «W /(«)** (27) 



m 2 \NJ 

m—1 



with 



2/ n , r , 12 



Cin(47rx) H sin(47ra;) — 1 



(28) 



.r- 

It is straightforward to check that f(x) = Oq(x 2 ) and f(x) = Oqc ( ^ n< ^2 1 so the integral in (127|) is perfectly 
well defined. 

For the second sum, the power N — 2 prevents us from using a similar technique. However, we can under- 
stand qualitatively why this term is subdominant. The factor s (S) quickly goes to zero at large TV because 
s(x) = sin(27rx)/(27Tx) < 1 except when m/N -c 1. Therefore, only the modes m -C TV contribue to the sum but 
as TV grows, the number of such modes also grows. However, becaues in the limit x — > the function ^ 2 (x)/x 2 
goes to zero, their contribution tends to zero and therefore the total sum goes as TV times a term of limit zero. 
A similar argument shows that the third term is also subdominant. 

As a conclusion, in the large TV limit, we have 

SH^(?f (29) 

The integration over k is trivial since is independent of k and yields 



f(u)du N (30) 



The integral can be computed numerically to obtain our final result for the power radiated by loops with 
TV ^> 1 left-moving and right-moving kink^l 

(P)^Gfj 2 T'N with r'«13 (31) 



3 If we consider loops with N left-moving kinks but only 2 right-moving kinks, a similar computation leads to 
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2.4 Numerical study 

2.4.1 Computation for a given loop 

We now explain how the power radiated by a given loop belonging to the class defined in Section 12.11 can be 
computed numerically. This calculation was performed analytically by Garfinkle and Vachaspati in [21j in the 
case N = 2. In this section, we generalize their method to an arbitrary value of N. The infinite sum over m in 
([7]) can still be performed analytically. The integral over k is done numerically over a grid. 
Let us start by rewriting (JTUJ) in the following way 



where we defined 



. N-X 



n=0 



(32) 



2tt 
~N 



±(n + 1) - k • A ±i . 



and 



2tt 

iV 



±n 



(33) 



as well as 



Inserting this into 
the unit vectors A± iTl 



D i A ±,n 

±>" k-A ± ,„Tl' 



(34) 



and l[7|). we can write the power radiated in the direction k as a quadruple sum over 



dP 4G/1 2 (i j) Q p ) 1 / im(-a+-a _ , +a++ a - ) . \\ ,„_, 

dk = ^T~ E A ij,ip D +,n D An' D +, q D -, q > I 2^ ^2 ( e " " 9 + similar terms J I . (35) 

n,n' ,q,q f — \m— 1 / 

All 16 terms in the parenthesis above come from the expansion of (e~ ima ™ _ e - imo 4)(e~ lmQ »' -e~ lma ™')(e lma ^ - 
e "«"J"jj e "" a ,' _ e mia ,,'^ so they are all of the form ±e"" 7 . Furthermore, ^ has to be real so we can replace all 
terms like e lmi by cos(m7). We can now use the following summation formula to simplify ([33]) 

V> cos(m 7 ) 1 2 7T 2 

V 2^ — = - ( 7 mod 2tt -7T - — . 36 

' m z 4 12 

m— 1 

Since 8 of the 16 terms have a plus sign while the other 8 have a minus sign in front, the contributions from 
— j2 in (|36p cancel in (|35|) and we are left with 



dk ~ TT 3 ^ 1X n,lv U +,n U -,n' U +,q U -,q' 

n,n' ,q,q' — 



2 

((— a^— a~,+a£+a~,) mod 2-k)— tt J +similar terms 



(37) 



In the cases where the direction of emission is exactly equal to some A +j „ (resp. to some A_ i<?0 ), this 
expression has to be corrected by removing all terms with n = no or n' = no (resp. q = qa or q' = q$) in the 
quadruple sum. 

The number of operations that one has to perform in order to compute numerically 4? for any given 
direction k grows as 0{N A ) and quickly becomes prohibitive. The computation can be made faster by using 
some symmetries but N of the order of 100 seems a reasonable limit if one wishes to perform the integration in 
([T2"j) to obtain the total power radiatcc(f|. 



Y^Zi 7^ [ a + {~l) m s N/2 - V] Cin(27rm) instead of Eq. J25} which can be shown to grow only logarithmically with TV. 
This is consistent with the results of |26| that small scale structure only traveling in one direction on the string does not lead to 
efficient energy loss. 

4 As an alternative to our procedure, we could have used the results of 1231 with the particular advantage of getting rid of the 
numerical integration, both a potential source of error and a time consuming step. We only chose to stick to our direct generalization 
of 12 1 1 because it is simpler and because the smoothness of the power as a function of k enabled us to perform the integral with a 
reasonably small grid. 
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2.4.2 Results 



Combining the procedure described in Sec. 12.3.11 to draw random loops in our class and the method of the 
previous paragraph to compute their power, we can now statistically estimate (P) for various values of N. The 
results, shown in diamonds in Fig[T]and well fitted by the solid (red) line, confirm the linear behavior derived in 
Sec. 12.31 when the number of kinks becomes large. The best fit of our data with a linear function yields a value 
for the slope equal to 13.5, showing good agreement with (|3ip. We also observed that the ratio of the standard 
deviation to the average value decreases as N grows so that the average value should be a good estimator of 
the typical power radiated. 




Figure 1: Average power radiated by a loop with N kinks in units of G/i 2 as a function of N. The diamonds 
correspond to the distribution of loops defined in Sec. I2.3.1l for which the distribution of kink angles is broadly 
peaked around 7r/2 and are very well fitted by the red straight line illustrating the linear growth at large N 
found in Sec. 12.31 The top and the bottom set of points (and the associated dashed blue and dotdashed black 
straight lines) correspond to piecewise constant loops with fixed kink angles of respectively 135° and 7r/10, 
showing that the power radiated increases with the amplitude of the kinks. 

The distribution defined in Sec. 12. 3. H is flat in the sense that each of the A± )7l has a uniform probabilty on 
the unit sphere (or rigorously half of them, in order to obtain a closed loop) . This induces a whole distribution 
of kink angles, broadly peaked around the typical value n/2. 

We would now like to know how the amplitude of the kinks affects the power radiated. We therefore modify 
the procedure to restrict ourselves to loops with fixed kink angles in the following way. Given a and two random 
unit vectors A± ; o, we randomly draw each unit vector A±. n +i (for n between and N/2 — 2) with uniform 
probability among those that satisfy A.± tTl ■ A± i7l +i = cos a. With the usual definitions in the literature (see for 
example [U]), this corresponds to a kink amplitude (or sharpness) A = i||A± >n — A± j?1+ i|| = sin(a/2). Note 
that strictly speaking, not all of the 2N kinks on these loops have an angle a but only 27V — 4, the other 4 
being potentially significantly different. 

The results are shown in Fig[T] for two different values of a, respectively 135° for the upper set of points 
(boxes) and n/10 for the lower set (circles )@ Here again, we find that the average power increases linearly with 
the number of kinks. We note however that the standard deviation becomes important at small a so for a given 

5 The value 135° was chosen in order to compare our results with those of 1241 . In addition to their argument based on the 
quadrupole radiation formula, the authors performed numerical simulations of loops with 8, 16 and 32 kinks with opening angles of 
135° to compute the numerical proportionality constant (r^ in their notation). Note that these were simulations of the dynamics 
of Nambu-Goto loops including gravitational backreaction, an approach completely different from the numerical computations of 
the radiated power that we performed in this section. Taking into account the slightly different definition of the word kink used 
in their paper which introduces a factor of 2, their results translate into a value rj 25 — 50 for the slope corresponding to this 
particular value of the kink angle. 
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loop, the actual power radiated can be significantly larger than the average valued The best fits give slopes 
equal to 24 for a = 135° and 2 for a = 7r/10. As could be expected, we find that the power increases with the 
kink angle. 

2.5 Application to loops with junctions 

We now conclude this section with a few remarks on the applicability of our results in the context of kink 
proliferation as described in [18j which occurs on non-periodic loops containing junctions and gives rise to a 
whole distribution of kink angles. First, we argue that although we have worked here only with closed periodic 
loops, the periodicity of the functions x± does not play an important role and as a consequence, we expect our 
results to hold for the loops with junctions. Indeed, the fact that x± are periodic induces correlations within 
each set of vectors A +j „ and A_. n (in our class of loops, this translates into the choice A.„ t ± = — A n _jv/2.±) 
but these correlations only give rise to the second and the third term in (|26|l and those are subdominant at 
large N. 

A more problematic issue is that of the value of T' corresponding to the complicated distribution of kink 
angles produced during kink proliferation. The population of kinks resulting from proliferation can be divided 
into a number k' of large amplitude kinks which dominate the signal received by an observer from a whole 
network of loops and a huge number of unobservable very small kinks. The tools that we developed in this 
section do not enable us to compute V for such an intricate distribution: because it is only known numerically, 
no analytical approach like that of Sec. 12.31 is possible and the number of kinks exceeds by far the possibilities 
of the numerical approach of Sec. 12.41 However, we can still put a lower bound on the power radiated by such 
loops by considering only the large amplitude kinks. In this case, we have according to our results P ~ G^T'k 1 
with a value of V of the order of 10. Given that we will only be interested in order of magnitude estimates 
in the remainder of this paper and that several parameters already suffer from comparable uncertainties, we 
can replace V in the expression of P by the usual constant r ~ 50 expressing the power radiated by smooth 
loops in units of G/i 2 in order to simplify our expressions. This lower bound will be a reliable estimate of P 
if the contribution from the small amplitude kinks is negligible. If instead the power is actually dominated by 
the small amplitude kinks, then P will depend on their number rather than on k' . In this case, we can only 
conclude that P ^> TGfi 2 k'. Distinguishing between those two situations is not possible with the tools that we 
have developed so far and is a model-dependent issue since different types of strings may give rise to different 
kink amplitude distributions. As a consequence, in the rest of this paper, we will discuss the implications of 
both scenarios. The conclusions of this discussion can be summarized as follows: 

f P~ TG/i 2 • k' Power dominated by large kinks ,„.. 
I P> TGfj 2 ■ k' Power dominated by small kinks 

3 Consequences on the GW signatures from kink proliferation 

The effect of kink proliferation on the GW signal emitted by a network of cosmic strings was computed in |18j 
in the general framework of a network where two types of loops coexist: "standard" loops with a number of 
kinks of order unity, and loops with junction^] that have gone through a kink-proliferation phase. The average 
number of kinks per loop of the second subnetwork was taken to be a free parameter k' S> 1. 

Obviously, the predictions for the signal depend crucially on the number density of each type of loop since the 
rate of bursts is directly proportional to this quantity. In this section, we revisit the results of |18j , improving 
the estimation of the density of loops with k' kinks by taking into account the fact that they radiate more 
gravitational power than initially considered and therefore decay faster. In Sec. 13.11 we derive our new expression 
for the density of loops which is reduced at least by a factor k! . The consequences of this lower density on the 
gravitational wave bursts and on the stochastic background are studied in sections 13.21 and 13.31 respectively. 

6 This is due to the following: when a is small, the A_|_ n will tend to be gathered together, densely populating a small region 
instead of covering the whole unit sphere. Now, in the cases where A+ : o and —A—fi are drawn not too far from each other, these 
two regions intersect and the chances that one of the A_|_ iT1 becomes really close to one of the — A_ >n resulting in an extremely high 
instantaneous velocity and important radiation (an equality would result in a degenrate cusp as discussed in Sec. 12.21 will be high. 

7 In this whole section, slightly abusing language, we take this expression to include also the periodic daughter loops (with no 
junctions) that may be formed due to the possible unzipping of junctions. In particular, when we talk about the lifetime of loops 
with junctions, we mean the time of decay by gravitational radiation of the loops formed with junctions and of their possible 
daughter loops obtained by junction unzipping and not the time of unzipping itself. 
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3.1 Lifetime of loops with many kinks and density of loops 



3.1.1 One-scale model - standard case 

The predictions in [18] are based on a modified (in order to account for the presence of junctions) version of the 
standard one-scale model for the network evolution that was used for the original estimates El UHl H] . We start 
here by summarizing the key features of this model (see [27J for a detailed discussion, together with [5] for for 
instance the factors of p) . 

A string network is made of both long strings and subhorizon closed loops. The long string component can 
be characterized by the typical distance between long strings L st , the typical distance £ s * beyond which the 
directions of string segments become uncorrelated and the typical wavelength of the smallest wiggles l W i gg \ cs on 
the long strings. In the so-called scaling regime, all these lengths are proportional to the Hubble radius and 
we have L st ~ p 1 l' 2 t and £ s * ~ t. The value of biggies is the subject of intense debate in the community and is 
usually written in terms of the parameter a under the form Z W ig g ios ~ ext. Using a crude estimate of the efficiency 
of gravitational backreaction to round off the small scale structure on the long strings [25] gives a ~ TGfi where 
the constant is roughly equal to 50. 

At any given time, loops are chopped off from (self-) intersection of long strings. Those large loops then 
fragment into many small loops (which give the dominant contribution to the GW signal) and whose typical size 
I is that of the smallest wiggles on the long strings I ~ at. Under this assumption, the population of loops can 
be determined by remarking that the scaling regime requires a given rate of energy loss for the infinite string 
network (p -1 /ii per Hubble volume and per Hubble time) so the number of such loops formed in a Hubble time 
and within a Hubble volume has to be p fit/(fj,l) ~ pa. These loops then decay by gravitational radiation. 
The power radiated by these smooth loops has been computed in [Ml HOI [2T] and is proportional to G/i 2 with 
the same proportionality constant r ~ 50 as in the previous paragraph. The lifetime of loops formed at time t is 
therefore r ~ yg^ 2 ~ ^tg^ 2 ~ ^ n °th cr words, all the loops formed at a given time t decay in a Hubble time 
and with good approximation all loops have the same size at at any time t which leads to a distibution of loop 
lengths approximately equal to a delta function n(l, t) ~ -^^6(1— at) that we will write from now on n(t) ~ ^jt- 

Some authors have computed the gravitational wave predictions in different scenarios for the evolution of 
the loop network, for instance including smaller loops with a = eYG^i and e < 1 or large loops [H El HH]. For 
the sake of simplicity we will however stick, as in [18] , to a modification of this simplest scenario to analyze the 
impact of kink proliferation. 



3.1.2 One-scale model - with proliferation: old version 

Several authors have studied the evolution of a network of cosmic strings with junctions, both analytically 
and numerically [Ml EH EH [HI E31 EH EH ES] and found that the scaling behaviour of the infinite strings 
still exists in the presence of junctions. Based on these results, the one-scale model presented in the previous 
section was modified in [T5J to take into account the population of loops with junctions in the case of a network 
containing different types of strings of comparable tensions and reconnection probabilities. The simplifying 
picture adopted there in order to isolate the effect of kink proliferation was that despite the formation of 
junctions, all the quantities related to the infinite string network remained unchanged. 

We are interested in the population of small loops and in particular in the fact that now two types of loops 
can form: standard periodic loops and loops with junctions on them. Because both types of loops are produced 
as in the standard case by iterative fragmentation of larger loops, it was argued that their typical size was that 
of the smallest wiggles on the infinite network at. In order to describe these two different types of loops, a new 
parameter q defined as the fraction of loops formed with junctions was introduced. 

As in the standard case, loop production is the mechanism by which the infinite string network loses en- 
ergy and in order to ensure a scaling behavior, the total energy lost has to remain unchanged with respect to 
the standard case. Because both types of loops have the same energy, the total number of loops formed per 
Hubble time and within a Hubble volume remains \/{pa). However, this number now divides into (1 — q)/(pa) 
standard loops and q/ (pa) loops with junctions. Because the creation of a loop with junctions can only happen 
by fragmentation on a larger loop around a junction whereas the production of a standard loop can happen 
anywhere on the loop, we expect g < 1 and the exact value of this parameters is probably model dependent 
since it is related to the tendancy of a particular type of string to create junctions. When working with order of 
magnitude estimates, 1 — q can be replaced by 1 in the number of standard loops. Furthermore, in the scaling 
regime, q is expected to be time-independent. 
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We are interested in the number of kinks on each type of loop. The standard loops arc the product of the 
fragmentation of loops without junctions and therefore, as in the standard case, should have a number of kinks of 
order unity. On the second type of loops, the presence of junctions leads to a quick proliferation of the number of 
kinks and we call k! the average number of large amplitude kinks on these loops. Note that, as discussed in more 
detail in SectionHJ because the proliferation of kinks is exponential, k! is roughly equal to the number of kinks at 
the end of proliferation regardless of how the the duration of proliferation compares with the lifetime of the loop. 



The last step to obtain the number density for each type of loop is to take into account the decay of these 
loops by gravitational radiation. At that point, [18j made the conservative assumption that the power radiated 
by both types of loops was the same as in the standard scenario (P = TGp 2 ) and in particualar that the power 
radiated by the second type of loops was not affected by kink proliferation. As in the standard scenario, this 
leads to a number density 

n(t)r 

for standard loops and 
for loops with junctions. 



pat 3 



pat 3 



(39) 
(40) 



3.1.3 One-scale model - with proliferation: revisited estimates 

Since the loops with junctions are created with an invariant length I ~ at ~ TG/it, then their initial energy is 
\iat = TGp 2 t. Combining this with our expression in (|38|) for the power radiated by those loops, we find that 
for models where the power is dominated by the large kinks, the loops decajd in a time 

TGfiH t 

Tl ~ TG^ ~ ~k r (41) 

In the case where the small kinks provide the main contribution to the loss of energy from the loop, this value 
becomes an upper bound according to (|38|) . In both cases, this is much shorter than the value r; = t that was 
used in [18] . Therefore, at any given time, only a fraction of all the loops considered in |18j will be present in 
the network and the number density of loops with many kinks is at most 

k pat^ 

We now analyze the effect of this smaller density on our results. 



3.2 Rates and individual bursts 

The goal of this subsection is not to compute in detail the rate of bursts but rather to explain qualitatively 
the consequences of the shorter lifetime on the predictions for the GW signal. Therefore, instead of explicitely 
describing the dependence of the rate of kink bursts on the various parameters and on redshift (see sections 3 
and 4 of [TS] for more detail), we will only be interested here in the intuitive fact that it is proportional to the 
density of loops and to the number of kinks per loop. Since all the other factors entering the rate are equal for 
both types of loops, this approach will be sufficient for our purpose of comparing the contribution from standard 
loops and from loops with junctions in our model. 



Following the results of 13.1.21 it was found in [TH] that 

' functions \ q ■ k' ■ p -1 a -1 t -3 



/ A*j unc tj ons \ 

V A s tandard / old 



p 1 a 1 t 



qk' (43) 



with the conclusion that for models satisfying the condition qk' 3> 1 , the main contribution to the rate of kink 
bursts came from the loops with junctions. 



8 This simple calculation implicitely assumes that the number of kinks remains constant throughout the evolution of the loop. 
Though this is not necessarily the case, the analysis of Sec. [4] shows that because the number of kinks grows exponentially, our 
estimate for the lifetime of the loop is valid. 
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Now using our revisited one-scale model taking into account the faster decay of loops with many kinks, we 
find that this fraction is smaller by at least a factor k' 



f unctions \ < q-U-p-l a -H-*W) ^ 

V ^standard / new p 1 a 1 t 3 

and because q < 1, the contribution from standard loops always dominates the rate of kink events, making imme- 
diately clear that the original predictions made by 119| concerning the observability of individual bursts apply. 

Note that this leaves the final conclusions of [IB] on the observability of individual kink bursts unchanged. 
Indeed, it was argued there that the usual procedure to compute the rates and amplitudes had to be revised 
by introducing a lower redshift cutoff z c below which no statistical approach was valid due to the lack of loops 
in this regiorQ This led to the conclusion that the presence of kinky loops could only yield amplitudes for 
kinks bursts larger than the standard predictions in [19] in specific regions of parameter space where pulsar 
constraints already forbid this anyway. 



3.3 Stochastic background 

For any given type of burst, the rate of events coming from high redshifts is so large that an observer cannot dis- 
tinguish the individual bursts and only sees their incoherent superposition as a stochastic background. Damour 
and Vilcnkin [19j showed that the sum over the overlapping bursts could be written in terms of a redshift 
integral of the square amplitude of individual bursts emitted at redshift z weighted by the number of bursts 
from that redshift. The crucial point is that one has to start the integral at the redshift Zb^b that corresponds 
to the minimal rate for the bursts to actually overlap at the detector. The integral ends at Zhf beyond which 
the loops are so small that they only emit at frequencies larger than the frequency of observation. Obviously, 
this calculation is only valid for redshifts larger than z c . This reads 

hl(!) = h 2 (z)f- 1 dN(z). (45) 

Our goal here is to investigate the modifications with respect to the results of [TS] induced by the smaller 
density of loops with junctions n' . Three of the quantities appearing in (|45[) depend on n', namely z c , £&->6 
and dN(z) (the rate of bursts coming from the redshift interval [z, z + dz]). However, both in the case of the 
kink background and of the kink-kink background (superposition of bursts emitted at kink-kink encounters on 
the strings), the integral is always dominated by the largest redshift Zhf (see Appendix B of [IB]) which only 
depends on the size of the loops and therefore remains unchanged. As a result, the only corrections come from 
the factor dN(z) which is directly proportional to n' in the case of loops with junctions and now smaller by 
a factor k! or larger. This means that the final formulas for the amplitude of the kink and of the kink-kink 
background (respectively formulas (53) and (57) of [IB]) need to be divided by at least a factor \/k' . 

Let us start by discussing the kink background. We now have, for the contribution from loops with junctions 
in the case where the large kinks dominate the power that they radiate, 

(\ !/3 
1 + J , (46) 

the right hand side being an upper bound for /i kmk in the case of a small kink dominated power. The contri- 
bution of the standard loops is obtained by substituting q = 1 in the previous equation. Since q < 1, we now 
find that the kink stochastic background is actually always dominated by the standard loops, even in the models 
where the product qk' is larger than 1. 

We now turn to the kink-kink background for which our conclusions will be different depending on whether 
the power radiated by the loops with junctions is dominated by the large kinks or by the small ones. 

In the first case, the contribution from loops with junctions obtained in [18j has to be divided by \[k! and 
now reads 

1 /2 

h^ k ~ {qk'f /2 ■ 10(G^)aVV 1/2 e (f - -i- J (1 + «*o/r 3/2 ( 1 + ^) ■ (47) 

\ at J \ z eq J 

9 Also note, though this is not necessary here to conclude that the dominant kink signal comes from the standard loops, that the 
smaller redshift cutoff applying to loops with junctions is now larger than in |18| by a factor (fc') 1 / 3 or larger. 
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Here, the contribution from standard loops is recovered by setting q = k' = 1 and so, the conclusion that it is 
dominated by loops with junctions provided that qk' > 1 remains valid. Furthermore, it is easy to check that in 
this case the k-k bacground remains larger than both the cusp and the kink backgrounds computed in the stan- 
dard case. Indeed, we now have h k k ( - ncw - ) ^ h k k ^ old ' /y/W and h k ( - ncw - ) n*> h\ ^ oXA ' /\fqW. On the other hand, it 
was shown in [TH] that he k /h k ( ' old ' ) ~ \fk' '. Combining these equations leads to h k k ^ ncw - ) //ic ^ new - 1 ~ y 'qk' '. 
Note that h k < - now < 1 corresponds to the amplitude of the standard kink background, and that, according to [55] , 
it is of the same order of magnitude as the standard cusp background. 

The case where the small kinks dominate is more complex because the estimate in (|47[) now becomes an 
upper bound and it is impossible to compare the contributions from standard loops and loops with junctions as 
well as the backgrounds from different types of bursts. For now we simply notice that if the radiation from the 
small kinks is too large, then even in models with qk' > 1, the kink-kink background will not be the dominant one. 

From an observational point of view, the conclusions of this section arc that the amplitude of the kink-kink 
background is reduced with respect to the results of [TS] by a factor y/W at least and that the region of parameter 
space for which the dominant background remains the kink-kink could be smaller than computed in [18] . 



4 Kink proliferation versus power radiated: towards estimating k' 

Our results on the power radiated by a loop with many kinks enables us to discuss the order of magnitude of 
the parameter kl entering the observational predictions of Sec. [3J 



4.1 Unspoiled exponential kink proliferation 

Let us consider a loop with junctions of typical length L undergoing kink proliferation due to the transmissions 
and reflexions of kinks propagating through the junctions. We neglect for now all the potential mechanisms 
identified in [T7] that could alter the course of proliferation, namely the self-intersections of the loop which could 
chopp off smaller loops without junctions but with many kinks, the unzipping of junctions and the gravitational 
backreaction (by this we mean the deviations to the Nambu Goto dynamics on a fixed flat background and in 
particular the rounding off of the kinks; as explained below, we however take into account the fact that the 
energy of the loop is indeed decreasing). We discuss the effect of these additional physical ingredients in Sec. 
14.21 Under these assumptions, it was shown in [17j that the number of sharp kinks N(t) on the loop grows 
exponentially: 

JV(t)~ P t/L (48) 

where the origin of time was set at the formation of the loop. The coefficient p > 1 depends on the tensions of 
the strings in the loop. In the case where all the strings have the same tension, p is maximal and roughly equal 
to 2. 



As the number of kinks increases on the loop, it radiates more power until, at t = Tf it has lost all its initial 
energy which is roughly equal to pL. 

We start by discussing the case where the power radiated by the loop is dominated by the large amplitude 
kinks during the whole course of proliferation. Using (|38|) . we find that between t and t + dt, the loop looses an 
energy dE ~ TGp 2 ■ N(t)dt. Therefore, tj satisfies the equation 

,2 



TGp 2 ■ I N{t)dt ~ pL (49) 

which, combined with (j4"5]) leads to 



provided p is not too close to 1. 



How is this related to the parameter k' appearing in the computations of the gravitational wave signal in 
Sections 13.21 and [ 3.31 for instance? The number of kinks does not remain constant throughout the lifetime of the 
loop and k' has to be defined as the average number of kinks per loop which is equal to the time average of the 
number of kinks on a given loop. Therefore, k' will be somewhat smaller than N{jf) but, given the exponential 
growth of the number of kinks, we expect it to be comparable. This gives us our final estimate for k' 

k' ~ (TGp)- 1 . (51) 
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In the case where the small kinks actually provide the main contribution to the loss of energy from the loop 
for at least part of the proliferation phase, our estimate turns into an upper bound 

This is the number of kinks one obtains by assuming that proliferation occurs in the ideal conditions men- 
tioned earlier. In other words, when the only limitation to proliferation is the decay of the loop by gravitational 
radiation. In this case, the duration of proliferation is equal to the lifetime of the loop. In |17j . various mecha- 
nisms which might be responsible for an alteration of the proliferation of kinks were proposed. We now discuss 
qualitatively how these can affect our estimate of k'. 

4.2 Discussion of other effects 
4.2.1 Unzipping of junctions 

The first mechanism is related to the potential unzipping of the junctions. In particular, simulations of realistic 
loops with junctions in |17| showed that one of the segemcnts of string typically ends up shrinking to a point 
and at that time, both junctions collide. What happens at that point in not encoded in the effective modified 
Nambu action that was used and certainly depends on the underlying theory 37, 38J but one possibility is that 
the junctions disappear or unzip to release two periodic closed loops. In that case, proliferation brutally stops 
and the number of kinks then remains constant. It is interesting to remark the typical time for this collision 
to happen (a few L according to our simulations) is comparable with (and seems to be slightly smaller than) 
the lifetime predicted by considering the gravitational decay which, as can easily be seen from (|48p. is given by 
In (TGn)^ 1 in units of L. Therefore, the number of kinks can be significantly smaller than that predicted in (|51[) . 



4.2.2 Self-intersections 

Another physical phenomenon that we neglected is the possibility for the segments of string to self-intersect 
during the evolution and reconnect to chopp off small loops carrying many kinks which therefore do not enter the 
proliferation anymore. This effect has not been studied quantitatively for loops with junctions. In particular, 
it is not clear to what extent the frequency of these self intersections will increase as the number of kinks 
grows. Moreover, the effect will depend on the value of the reconncction probability since a low value of p 
would suppress the production of small loops even when many self intersections take place. In any case, as 
we explain below, we do not expect that multiple rcconnections affect drastically the estimate for kl . Indeed, 
when a loop without junctions is chopped off, it withdraws kinks from the loop where proliferation takes place 
but also instantaneously reduces its length. This has two effects. It increases the rate of proliferation as the 
exponent in (|48p is now larger but it also decreases the energy of the loop and therefore its lifetime and the 
duration of proliferation. 

Let us consider a very simplistic toy model to measure the importance of these competing effects. We start 
with a loop containing three segments of length L meeting at two junctions. Suppose that at some instant to, 
each of the three segments self-intersects and has its length reduced by a factor a < 1 (the length of the loop 
released is then (1 — a)L). Assuming that the kinks are evenly spread over the segment, the number of kinks 
on each segment right after reconnection is Nq = a p*°/ L and from then on, it evolves as N(t) = NQp^ t ~ t °' l ^ aL \ 
Solving for the number of kinks at time Tf when the loop has lost all its energy by gravitational radiation, we 
recover the same result as in the case with no reconnection N(rf) w (TGfi)^ 1 + No ~ (TGfi)^ 1 . The value of 
k' is slightly increased by the presence of kinks in the three loops produced at the reconnections but as there is 
no proliferation on these, we expect this contribution to be subdominant. 



4.2.3 Backreaction 

While we accounted for the loss of energy by gravitational radiation in our analysis, we completely neglected 
gravitational backreaction, namely the fact that this loss of energy causes the dynamics to deviate from the 
Nambu action on a fixed background. In our context, backreaction has at least two important effects. The first 
one is to round off the kinks SO] , with a timescale that is probably of the order of the lifetime of the loop 
Tf that we computed by looking at radiation. This tends to reduce the number of kinks in the system but at 
the same time, it increases the lifetime of the loop and therefore the duration of proliferation. We find however 
unlikely that this will contribute to increase the number of kinks with respect to (|51|) . The second effect is that 
the string shrinks as it loses its energy, or in other words, its length is actually a (decreasing) function of time. 
This will certainly alter the p l l L behavior that was derived in [T7] where L was the initial length of the segments 
of string, in particular towards the end of the evolution where the rate of proliferation will become increasingly 
large. 
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5 Conclusion 



In this paper, we have investigated the loss of energy by gravitational radiation from cosmic string loops with 
many kinks. Using analytical tools, we found that the power radiated by piecewise linear loops with N kinks and 
a distribution of kink angles roughly peaked around the value tt/2 is proportional to N with a proportionality 
constant of the form T'G/j 2 and the order of magnitude V ~ 10. We also investigated how our result depends 
on the typical amplitude of the kinks by computing numerically the power radiated by loops with fixed kink 
angles. We showed that the linear behavior in the number of kinks is robust and that V is an increasing function 
of the typical kink angle. This is the first time that radiation from loops with many sharp kinks was studied 
using a full relativistic formalism. 

We then applied our results in the context of kink proliferation on loops with junctions, a phenomenon de- 
scribed in [17j . The distribution of kinks on loops with junctions is complex and hard to study even numerically, 
in particular because it contains a huge a number of very small kinks. For this reason, precisely estimating the 
power radiated by such loops is impossible. However, our results from Sec. [2] enabled us to determine a lower 
bound on this power by considering only the interaction of the k' sharp kinks that dominate the contribution to 
the gravitational wave burst signal received by an observer from a whole cosmological network of loops. More 
precisely, we found that for loops with junctions, P > TGfj, 2 ■ k' where r ~ 50 is the standard value of the power 
radiated by smooth loops in units of G/x 2 . For models where the contribution from the small kinks to the power 
radiated is subdominant, this lower bound will be saturated. In any case, the power radiated by loops with 
junctions is much larger than the estimate P ~ YG^l 2 used in [TS] . 

As explained in Sec. this affects the observational predictions for the gravitational waves emitted by a 
network of cosmic string loops with junctions: since they radiate more, those loops decay faster and so their 
number density is smaller thereby decreasing the rates of events. We found that the rate of individual kink bursts 
should not be larger in networks with junctions than in standard networks. This strengthens the conclusion of 
[18] that, contrary to what one might expect at first sight, the presence of loops with junctions in a network and 
the associated kink proliferation do not provide any distinctive signature as far as the observation of individual 
bursts is concerned. 

According to [18j . the main observational consequence of kink proliferation is the existence of a strong 
stochastic background of gravitational waves resulting from the superposition of the large number of bursts 
produced at kink-kink encounters. More precisely, it was shown that for models where the fraction of loops with 
junctions q and the number of sharp kinks on them kl satisfy the condition qk' 3> 1, this kink- kink background 
dominates over the cusp and the kink one. Our analysis slightly weakens this conclusion: the characteristic 
amplitude of the kink-kink background is suppressed by at least a factor Vk' and so the region of parameter 
space in which it remains dominant may be reduced. We note however that for models where the total radiated 
power from the loops with junctions is dominated by the sharp kinks, the condition for the kink-kink background 
to dominate remains qk' 3> 1. 

Finally, our result enabled us to put an upper bound to the parameter k' resulting from kink proliferation 
by computing the decay of a loop by gravitational radiation as the number of kinks grows. By assuming that 
exponential proliferation lasts until the loop has radiated all its initial energy away, we found that k' < (TGfi) . 
We also discussed qualitatively how other effects could affect this upper bound but ultimately, a quantitative 
analysis combining all these effects would be necessary to provide a complete picture of the physics of kink 
proliferation and to estimate k', a crucial value for the observational predictions. 
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A Expression of 




in terms of s, \i and a 



In this appendix, we detail the calculation leading to Eq. (|23p . Our starting point is formula (|10[) 




JV-l 



(52) 



where we recall that 




(53) 
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This yields 

AT-1JV-1 

(w) = ^EE *Hn> (54) 

n=0 n'=0 

with the definition 

S^ n , = /a^e'^K^ As ~^°=<> As ) \ e ^(«-«'). (55) 

Depending on the values of n and n', some of (or all) the vectors appearing inside the average brackets are 
independent so that the average conveniently factorizes. In order to make the computation simpler, and to 
reduce the number of different situations that we need to distinguish, we will use the symmetry property 

Sn>n = f^n're'J to re duce the sum to terms with n < n' . This leads to 

AT-l 

™ 2 = E + E + ( 56 ) 

n=0 0<n<n'<JV-l 

The terms in the first sum simply reduce to 

S% n = Kc£). (57) 
The second sum only contains terms with n' > n, which we can rewrite 

S% n>>n = (ciCil lle^Ae^-^. (58) 
We need to distinguish the following cases: 



• n 1 = n + N/2: There are N/2 such terms. In this case, A„ = — A„/. All the factors between the brackets 



are uncorrelated except for C % n , C 3 *, and e *» m kA " . We then have 



n-l+N/2 

S^ +N/ , = i-l) m (ciK\ N /^^) II (^ k ' As ) (59) 

s—n+l 

• n' < n + N/2: All the factors between the brackets are uncorrelated except for C % n and g^Tp 1 ^™. i n this 

case, 

n'-l 

S^, <n+JV/2 = (^e^ k - A «} (<#) J] ( e ^^) e ^(«-«) ( 6 o) 

s— n+l 

• n' > n + N/2: Some of the terms under the product sign cancel out and all the factors that are left are 
uncorrelated. In this case, 

n—n'+N 

S:U>>n +N /2 = K)( C i>) II (e 2 ^^)^"-"') (61) 

s=n+l 

This leads us to define the following average values involving only one of the random vectors A n : 

P l = K) (63) 
P l = (Qe^ k ' A ") (64) 



<r« = (& n Ci*) (65) 

/ n K - , n +N/2 t 



^ = (Cici*e^^\ (66) 
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Given our procedure to draw the vectors A n , it is clear that none of the average values defined here actually 
depends on the index n. Note however that all of them depend on k and m but we keep this dependence implicit 
here. 

The expression of ) in terms of //•? , cr y , p\ p l and s reads 



m 2 



(l*+P+) = Na ij + t-^N ^s N / 2 ' 1 + ^*{s*) N / 2 -^ + K lP y* + K 2 j?ft* + K* 2 p l p>* (67) 
with 

N-l N-l 

k\ = Yl E e ^ {n ~ n>) sN ~ n ' ( 68 ) 

n=0 7i'=n+N/2+l 
N-l n+N/2-1 

K 2 = J2 £ e ^(n-«') s n'-» (69) 
n— n'— n+1 

In principle, we now need to compute the full tensors //•?, cr 1 - 3 , p\ p l as well as the value of s, K\ and K 2 - 
However the number of components that we will really need can be drastically reduced by using the symmetries 
of our problem. As pointed out earlier, since there is no preferred direction in our class of loops, (4^ ) will be 

independent of k and we had set k = (1,0,0). Now remember that we will eventually contract with 

the projection operator Aij i p which projects away any term proportional to k. Therefore, any component of 

(^I^I^*^ bearing an index 1 will be projected away. This means that we only need to calculate p 2 , p 3 ,p 2 , p 3 , 

a 22 , cr 33 , cr 23 , er 32 , p 22 , /i 33 , p 23 and p 32 . We can even go a little further by remarking that p z and p l can only 
be proportional to k z (the only vector in the game) so that 

p 2 = p 3 = p 2 = p 3 - (70) 

and the last three terms in (|67p will be projected away. This means that we don't need to compute K\ and K 2 . 
Furthermore, directions 2 and 3 are equivalent so we have 

a 22 = a 33 ee a, p 22 = p 33 = p, (71) 

a 23 = a 32 = p 23 = p 32 = 0. (72) 

which only leaves us with the three quantities a and p and s to compute. For the sake of clarity, we postpone 
their computations to Appendix IB1 

We can finally rewrite (|67|) in terms of s, p and a to recover Eq. (|23p 

( 7 !4*} = ^2 ( CT + (-^) m s N/2 - 1 p) - 5 a ) + Q y ' (73) 
where the tensor Q 1 ^ will be projected away by Aiji p . 



B Computation of s, fi and a as functions of m and N 

In this appendix, we detail the calculation of the functions of s, p and a defined in (|62j) and (|TTj) . Each of 
them is an average value of a quantity that depends on only one of the unit vectors drawn randomly on the 
unit sphere. We parametrize this unit vector using spherical coordinates (0,</)), so that the vector components 
are (cos(#), sin(#) cos(</>), sin(#) sin(^)). Since we assumed a uniform distribution over the unit sphere, <f> has to 
be uniformly distributed in the interval [0, 27r] while cos# is uniformly distributed in the interval [—1,1]. The 
average values now rewrite as double integrals over 9 and <p. 
Let us start with a 
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where we have used the variable change u = ^(1 — cos 0) to compute the integral over 9. In the last line, we 
have used the special cosine integral defined as 

Cin(x) - P 1 ~° 0St dt. (78) 
Jo t 



In the same fashion, we can compute 

2* 

i+N/2 



V = (ClCl\ N/2 e^^) (79) 



(k- A„-l)(k- A„ + l) 



l-k-A„) _ ^ ^(1+k.A.) _ ^ e ^k.A„ \ (8Q) 



— cos 2 <j>d(f> sin6-, 7—^7 t — - T sinf^(l-cos6') N )sinf^ : (l + cose)')^ (81) 

4?r 7 7 (cos6»- l)(cos6> + l) V JV v V V TV v ') y ' 

/2?rm\ JV {2mn\ , . 

cos sin . (82) 



Finally, we have 



2tt 



4-7T ./q ./n 27rm \ JV 



JV /2ttto\ 

Hit (83) 



Of course, one can check ([70]) . (j7Tj) and (j72")l using this procedure and variable changes adapted to the 
symmetry. 

References 



[i] 

[2] 
[3] 
M 
[5] 
[0] 



[9 
[10 

[11 



T. W. B. Kibble, "Topology of Cosmic Domains and Strings," J. Phys. A9 (1976) 1387-1398 



E. J. Copcland and T. Kibble, "Cosmic Strings and Superstrings," 



Proc.Roy.Soc.Lond. A466 (2010) 623-657 arXiv : 091 1 . 1345 [hep-th] 



T. Damour and A. Vilenkin, "Gravitational wave bursts from cosmic strings," 



Phys. Rev. Lett. 85 (2000) 3761-3764| |arXiv : gr-qcT/ 0004075 



X. Siemens et at, "Gravitational wave bursts from cosmic (super)strings: Quantitative analysis and 
constraints," \Phys. Rev. D73 (2006) 105001[ |arXiv : gr-qc/0603115| 

T. Damour and A. Vilenkin, "Gravitational radiation from cosmic (super)strings: Bursts, stochastic 
background, and observational windows," \Phys. Rev. D71 (2005) 063510[|arXiv:hep-th/0410222 



M. G. Jackson and X. Siemens, "Gravitational Wave Bursts from Cosmic Superstring Reconnections," 



JEEP 06 (2009) 089 arXiv: 0901 .0867 [hep-th] 



E. O'Callaghan, S. Chadburn, G. Gcshnizjani, R. Gregory, and I. Zavala, "Effect of extra dimensions on 
gravitational waves from cosmic strings," Phys. Rev. Lett. 105 (2010) 081602. IT003 .43951 

E. J. Copeland, T. W. B. Kibble, and D. A. Steer, "Collisions of strings with Y junctions," 
\Phys. Rev. Lett. 97 (2006) 021602] |arXiv:hep-th/060li53[ 

E. Copeland, T. Kibble, and D. A. Steer, "Constraints on string networks with junctions," 
\Phys.Rev. D75 (2007) 065024[ |arXiv:hep-th/06 11243 [hep-th] | 

E. Copeland, H. Firouzjahi, T. Kibble, and D. A. Steer, "On the Collision of Cosmic Superstrings," 



Phys. Rev. D77 (2008) 063521 arXiv : 0712 . 0808 [hep-th] 



N. Bevis and P. M. Saffin, "Cosmic string Y-junctions: A Comparison between field theoretic and 



Nambu-Goto dynamics," Phys.Rev. D78 (2008) 023503 arXiv: 0804 .0200 [hep-th] 



18 



[12] 

[13] 
[14] 
[15] 
[16] 
[17] 
[18] 
[19] 
[20] 
[21] 
[22] 
[23] 
[24] 
[25] 
[26] 
[27] 
[28] 
[29] 
[30] 
[31] 
[32] 
[33] 



P. Salmi, A. Achucarro, E. Copcland, T. Kibble, R. de Putter, et ai, "Kinematic constraints on formation 



of bound states of cosmic strings: Field theoretical approach," Phys.Rev. D77 (2008) 041701| 
|arXiv: 0712 .1204 [hep-th]| 



P. Saffin, "A Practical model for co smic (p,q) superstrings," JEEP 05 09 (2005) 011 
|arXiv:hep-th/0506138 [hep-th]| 



E. J. Copeland and P. Saffin, "On the evolution of cosmic-superstring networks," JEEP 0511 (2005) 023 
|arXiv:hep-th/0505110 [hep-th]| 

A. Rajantie, M. Sakellariadou, and H. Stoica, "Numerical experiments with p F- and q D-strings: the 



formation of (p,q) bound states," JCAP 0711 (2007) 021 arXiv : 0706 . 3662 [hep-th] 



P. Binetruy, A. Bohe, T. Hertog, and D. A. Steer, "Gravitational Wave Bursts from Cosmic Superstrings 
with Y- junctions," \Phys. Rev. D80 (2009) 123510[ |arXiv: 0907 .4522 [hep-th] [ 

P. Binetruy, A. Bohe, T. Hertog, and D. A. Steer, "Proliferation of sharp kinks on cosmic (super-)string 
loops with junctions," |arXiv: 10 05 . 242 6 [hep-th] | 

P. Binetruy, A. Bohe, T. Hertog, and D. Steer, "Gravitational wave signatures from kink proliferation on 



cosmic (super-) strings," Phys.Rev. D82 (2010) 126007 arXiv : 1009 . 2484 [hep-th] 



T. Damour and A. Vilenkin, "Gravitational wave bursts from cusps and kinks on cosmic strings," 



Phys. Rev. D64 (2001) 0640 08, arX iv:gr-qc/010 4026. 



T. Vachaspati and A. Vilenkin, "Gravitational radiation from cosmic strings," 
\Phys. Rev. D 31 (1985) no. 12, 3052-3058| 

D. Garfinkle and T. Vachasp ati, "RADIATION FROM KINKY, CUSPLESS COSMIC LOOPS," 
\Phys. Rev. D36 (1987)~2229} 

B. Allen, P. Casper, and A. Ottewill, "Analytic results for the gravitational radiation from a class of 



cosmic string loops," Phys.Rev. D50 (1994) 3703-3712 arXiv : gr-qc/9405037 [gr-qc] 



B. Allen and P. Casper, "A Closed form expression for the gravitational radiation rate from cosmic 



strings," Phys.Rev. D50 (1994) 2496-2518 arXiv : gr-qc/9405005 [gr-qc] 



J. M. Quashnock and T. Piran, "Effects of gravitational back reaction on small scale structure of cosmic 
strings," \Phys.Rev. D43 (1991) 3785-3788| 

S. Weinberg, Gravitation and Cosmology: Principles and Applications of the General Theory of Relativity. 
John Wiley & Sons, Inc., 1972. 

X. Siemens and K. D. Olum, "Gravitational radiation and the small scale structure of cosmic strings," 



Nucl.Phys. B611 (2001) 125-145 arXiv : gr-qc/0104085 [gr-qc] 



A. Vilenkin and E. P. S. Shcllard, Cosmic Strings and Other Topological Defects. Cambridge University 
Press, 2000. 

F. R. Bouchet and D. P. Bennett, "Millisecond pulsar constraint on cosmic strings," 



Phys. Rev. D 41 (1990) no. 2, 720-723| 



C. J. Burden, "Gravitational radiation from a particular class of cosmic strings," 
\Physics Letters B 164 (1985) no. 4-6, 277 - 281} 

S. H. H. Tye, I. Wasserman, and M. Wyman, "Scaling of multi-tension cosmic superstring networks," 



Phys. Rev. D71 (2005) 103508, arXiv : astro-ph/0503506. 



A. Avgoustidis and E. P. S. Shellard, "Velocity-Dependent Models for Non-Abclian/Entangled String 



Networks," \Phys. Rev. D78 (2008) 103 510 arXiv: 0705. 3395 [astro-ph] 



A. Avgoustidis and E. J. Copeland, "The effect of kinematic constraints on multi-tension string network 
evolution," \Phys.Rev. D81 (2010) 063517| |arXiv: 0912 .4004 [hep-ph]| 



M. Hindmarsh and P. M. Saffin, "Scaling in a SU(2)/Z(3) model of cosmic superstring networks," JEEP 
08 (2006) 066, |arXiv:hep-th/0605014[ 



19 



[34] J. Urrcstilla and A. Vilcnkin, "Evolution of cosmic superstring networks: a numerical simulation," 



JEEP 02 (2008) 037| |arXiv: 0712. 1146 [hep-th]| 



[35] M. Sakellariadou and H. Stoica, "Dynamics of f/d networks: the role of bound states," J CAP 0808 
(2008) 038. 10506 . 32T31 

[36] S. Olmez, V. Mandic, and X. Siemens, "Gravitational- Wave Stochastic Background from Kinks and 
Cusps on Cosmic Strings," \P~hys. Rev. D81 (2010) 104028[ |arX"iv: 1004. 0890 [astro-ph.CD]| 

[37] H. Firouzjahi, J. Karouby, S. Khosravi, and R. Brandenberger, "Zipping and Unzipping of Cosmic String 



Loops in Collision," Phys. Rev. D80 (2009) 083508 arXiv: 0907. 4986 [hep-th] 



N. Bevis et al, "Evolution and stability of cosmic string loops with Y- junctions," Phys. Rev. D80 (2009) 
125030, |arXiv : 0904 .2127 [hep-th] | 

J. M. Quas hnock and D. N. Spergel, "GRA VITATIONAL SELFINTERACTIONS OF COSMIC 
STRINGS," \Phys.Rev. D42 (1990) 2505-2520| 



[40] M. Hindmarsh, "Gravitational radiation from kinky infinite strings," Phys. Lett. B251 (1990) 28-33 



20 



